ID  ELEMENTS 


Simplest  type  of  FE  problems 
All  object  are  ID 
All  forces  are  ID 
All  stresses  /  strains  are  ID 


BASICS 


•  u=u(x)  :  Deformations 

•  8=8 (x)  :  Strain 

•  ct=o(x)  :  Stress 

•  f=f(x)  :  Body  forces 

•  t=T(x)  :  Tractive  forces 

•  P=P(x)  :  Point  loads 

•  a=Es,  e=du/dx 


Problem 


Discritization 


Element  connectivity  Matrix 
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:  Coordinates 
:  Local  coordinates 
:  Deformations 
(At  the  nodes) 


:  Deformation 
within  the  element 

Np  N2:  Shape  functions 


LINEAR  SHAPE  FUNCTIONS 


Nj  =(1-^/2 

N2=(l  +  ^)  /  2 


u  =  Nj  qj  +  N2  q2  =  (1  -  £)  /  2  *  (1  +  £)  /  2  *  q2 


Relationship  between  local  and  global  coordinates- 


4  =  -1+2*(x-x1)/(x2-x1) 


STRAIN: 


g  =  du/dx 

=  du/d£,  *  d^/dx 
u=N1q1+N2q2 

d£/dx  =  2/(x2-Xj) 
=  2/Le 

du/d£,  =  (-q1+q2)/2 
therefore- 


NjKl-^/2 

N2=(l+^)/2 


£  =  l/le  [-1  l][q,  q2]T 
=  Bq 


Where  B  is  the  element  strain  matrix  and 

B==  l/le  [-1  1] 

As  B  is  constant,  this  element  is  CONSTANT 
STRAIN  ELEMENT 

That  means  strain  inside  the  element  does  not 
vary. 

STRESS: 

a  =  EBq 


P.E.  approach 


n  =  X  [ a  T  e Adx  —  f  uT F Adx  -  f  uT Tdx  -  ^  uiPi 

L  L  L 


X  kJ  aT  eAdx  -sj  uT F Adx  -  11  uT  Tdx  -  E  Q,P. 


STRAIN  ENERGY 
FIELD 


FORCED  TRACTIVE  POINT 

FORCE  FORCE  LOADS 


Total  energy  of  the  body  is  the  sum  of  P.E,  of  all 
Elements. 


ELEMENT  STRAIN  ENERGY 


ke  is  element  stiffness  matrix 


CONTRIBUTION  OF  THE 
FORCES  TO  PE 


BODY  FORCES 


\qTNTfAdx  =  fAqT\NTdx 

J  Ntdx 
J  N2dx 


=  fAq 


1  Nxdx  -  area  under  the  triangle. 
=  (x2  - xx) / 2  =  Le  / 2 


T 

=  q  fe=  element  body  force  vector 


Thus,  body  force  on  element  gets  split  equally  at 
two  nodes. 


TRACTIVE  FORCES 


POTENTIAL  ENERGY 


Potential  energy  for  an  element  is- 

n.=XqTk.q-qTf.-qTT.-Y.qTP 

Total  Potential  energy- 

n  =zne  =  z  y2qTkeq-Y,qTFe 
F  _  \m.  |  Tie  fAie  n; 
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CONSIDER  A  SET  OF  ELEMENTS 


Element  No. 

Node  1 

Node  2 

1 

1 

2 

2 

2 

3 

3 

3 

4 

4 

4 

5 

5 

5 

6 

Q 


1 


Q 


6 


b,  ?Jb2  ?,rb,  «Jb4  qj 

ELEMENT  DISPLACEMENT  MATRIX 

Q=[q i  «2  ?3  -  -  -  9j 

GLOBAL  DISPLACEMENT  MATRIX 
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ELEMENT  STIFFNESS  MATRIX 


K= 


M 

h 

M 

h 

o 


o 


M 

A 


o 


,  ^4  ^4 


A  4 
^4 

4 

o 


4 

444._i_.A44 


o 


o 


i 


L 


2  *3 

44 


_44 

4 

44  44 


i 


l 


T 

_ 1 

1 

"  1 

^1 

- ► n  2 

^2 

T  3 

- ^  “ 

i 

Lr 

i+i 

GLOBAL  STIFFNESS  MATRIX 

1  2  3  4  5  6  .  n 
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where 


Vi  Vi  Vi 
K  IV  K  12’  k2V 

number  i. 


k'22:  stiffness  matrix  elements  of  element 


PROPERTIES  OF  K 


Symmetric 

Banded 

Can  be  sparse  (if  numbering  is  not 
appropriate) 

Is  N  X  N  (where  N  is  the  number  of  nodes 
in  a  1  D  problem) 


SPARSE  ‘K’  MATRIX 


THE  RESULTING  K  MATRIX:  SPARSE  AND  NON- 
BANDED 
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ASSEMBLING  GLOBAL  FORCE  MATRICES 
FROM  ELEMENT  MATRICES 


F  =  [(Fll  +Pj)(Fl2  +F12  +P2)(F22  +F31 +P3) 
(F32+F«+p4)...r 

^(f.+ro+p^F'.+p 

Where 


Feu  :  force  on  node  1  of  element  i 
Fe  :  force  on  node  2  of  element  I 

i2 

p.  :  point  load  on  node  number  i 


M+Ik  Ml+IL 
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ELEMENT  FORCE  MATRIX 
OR 

La  a] 

F=\fn  fn+fn  /22+/31  A  +  A 


GLOBAL  FORCE  MATRIX 


TOTAL  PE  USING  GLOBAL  MATRICES 

n  =  £l/2  *qT  *ke*q~YjqT  *  fe  -YsV7  *Te -HPi*Ui 

=YJll2*<iT*ke*q-Tj(iT*Fe-Tjpi*ui 

where 

Q  =  [qi  q2  q3  q*  ••••  qJT 

K  :  Global  stiffness  matrix 


F  :  Global  Force  matrix 


GALERKIN  'S  APPROACH 
| <yT  e  (</> )dV  -J aT fdV  - J t/>TdS  -  £  </>TP  =  0 

V  V  5 


<i>  (x ) 

G0  = 


J  (j) 

dx 


(f)  =  TV'F 
(^  )  =  5  ^ 

'p  =  i>,  *2r 


u 


Nq 

Bq 

EBq 


INTERNAL  WORK  DONE 


\<7t  s  (<p)dV  =  \qT BT EB'VAJx 


=  qT EA eBr B'V  J dx 

=  qT[EAJeBTB] ¥ 


qTkV 


Virtual  work  due  [  f  fdV  =  \'VT NT fA  dx 
To  body  forces-  *  ^ 

=  x¥TfAe$NTdx 

=  'F7/ ALj 

2  N  2d%  ~ 

_  *F  T  f  ^ 

J  2  L1. 

=  VTFe 


SIMILARLY  TRACTIVE  WORK  DONE 
\<pTTdx  =W  7 Te  where  Te  =  ^[l  l]r 

v  2 

Thus  Galerkin's  Equation  Becomes 

\crT  g  {(/>)dv  -  J  f  fdv  -  j <j)Ttds  -  ZfP  =  0 

v  S  S 

YJqTkeW-YJWTfe-YJx*TTe-WTP  =  0 
^TqTkeW-^TFe=0 

e  e 

I'l’  . '1\.. . '1' ,  | 


In  global  form- 


...  Galerkin '  s  approach 
...P.E.  approach 


TO  SOLVE  THE  PROBLEM 


Find  Deformations  Q.’s 


1.  Define  Boundary  conditions 

2.  Apply  minimization  of  PE 


Find  strains  s  =  B  *  q 
Find  stresses  a=Es=EBq 


SOLVING  FOR  O 


Ql^2\Q\  +  02^2202  +  02^2303  + . Ql^2nQn  + 


J 


qtkq-qtf 


BOUNDARY  CONDITIONS 

QPi  =  au  QP2  =  ai 

[e.g.  Q2  =0,  Q5=  0,  ...] 


-SINGLE  POINT  CONSTRAINT 


LET  US  CONSIDER  AS  AN  EXAMPLE 

Qi =  ai 

A  BECOMES 

n  =  ^  (a\K\  A  +  2Qo  +  a\K{  3C3  +  ...+a[K[nQn  + 

+Q2K2lal  +Q2K22Q2  +  Q2K23Q3  +  ...+Q2K2nQn  + 


QAa  +Q„Kn2a2  + . •+C„^sJ„2L) 

_0V*1  +  Q>A  +  03^3  +— +6,^,) 


MINIMUM  P .E .  PRINICIPLE 


—  (Klxax  +Q2K22  +  K23Q3  + . K2nQn  + 

Knax  +  Q2K22  +  K32Q3  + Kn2Qn )  -  F2  =  0 


—  (K3iax  +Q2K32  +K33Q3  + K3nQn  + 

^-na\  02^23  -^3303  . K3nQn  )  —  -^3  =  ® 


dn  =Q  ~(K;  1«1  +Q2Ki2  +Ki3Q3  + . +KinQn  + 

°Ql  La,  +Q2K2i  +  K3iQ3  +....KmQn)-Fl  = 


0 


Q2K22  +  Qt.K-2?,  + 

02^32  +  03^33  + 


+  QnK  2n  —  F2  Kllal 

+  Qn^in  =  Ft,  —  -^31^1 


Q2Ki2  +  Q3Ki3  + . +  QnKin  -  Ft  -  Kaax 
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K’=  OBTAINED  FROM  K  BY  DELETING  1st  ROW  AND  COLUMN 
Q’=  OBTAINED  FROM  Q  BY  DELETING  Q, 

F’  =  OBTAINED  FROM  F  BY  DELETING  F,  AND  SUBTRACTING  Kila1 
FROM  F, 


IF  INSTEAD  OF  Q,  =  a,,  WE  HAD  Q=a,  THE  SAME  STEPS  SHALL  BE 
CARRIED  PUT  BY  DOING  THESE  OPERATIONS  ON  ith  ROW  AND  COLUMN. 

REACTION  AT  THE  SUPPORT  (NODE  1  IF  Q,  =  a,) 

ELIMINATION  METHOD 


1 Q  +  ^1202  +  ^1303  +  *  *  •  *  +  KlnQn  -F1+R 1 
=KnQx  +KnQ2  +KuQ3  +....+KlnQn 


ENALTY  APPROACH 


ATTACH  A  SPRING  OF  STIFFNESS  C 

DEFLECT  THE  FIXED  END  BY  a, 
al 

DEFLECTION  OF  NODE  1-Q1 


wall 


AS  C  — >  oo,  Qx  — > 

-DEFLECTION  OF  SPRING  S  =  QX-  ax 


u,  =  y2  C(a  -  a,  )2  [i>.£.  0/  Sra/MJ] 


Q2K2]Q}  +  Q2K22Q2  +  Q2K22Q3  + ...  +  Q2K2nQn  + 


Ql^2\Q\  +  Q2K22Q2  +  Ql^liQi  +  •••  +  Ql^lnQn) 

~(Q\F\  +  Q2F  +  03^3  +  — +  0,A) 
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Multi  point  constraints- 

P\QP\  +  PiQpi  =  Ao 

for  example  -  ax=  2  a2 

we  can  solve  it  by  penaulty  approach 

nm  =  y2  qtkq + y2  c(Mpi  +  P2QP2  -p,)-qtf 

Since  C  has  very  large  value,  P.E.  takes  minimum  value 
When  (PjQ  j+  PjQpi-Po)  is  minimum. 


The  modified  stiffness  and  force  matrices  are- 


Rpl  pY^ pi  p2 

kp\p\  +  Cftl 

kp2p\kp2p2 

_kp2p\  +  CP\@2 
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Reactions  at  support  are  given  by- 

Rpi  -  ~C PiiPiQ pi  +  P2Q Pi  ~  Po) 
RP2  =  P i(P\Q P\  +  PiQ P2  ~  fio) 


STEPS  INVOLVED  IN  SOLVING  A 

ID  FE  PROBLEM 

•  Make  the  Geometric  Model 

•  Make  an  FE  Mesh 

•  Define  the  Loading 

•  Develop  Element  and  Global  Matrices 

•  Define  Boundary  Conditions 

•  Develop  Modified  Global  Matrices 

•  Solving  Using  Numerical  Techniques 


A  SIMPLE  PROBLEM 


P=1001b 


- 5725” — 

li=12” 

Elem  1:1-2 


E=3xl07  psi  12=12’ 

p=0.28361b/in3 

t=l”  Elem  2:2-3 


37T522 


k1  =  AE  /L  Tl  -1  1 

L-llJ 

k2  =  AE  /L  Tl  -1  1 

L-llJ 


f1  =  p  A  L  /  2  [1  1]T 


=  15.75  *  107/  12  |"l  -ll 

L-iiJ 

=  11.25  *  107/ 12  |"l  -ll 

L-iiJ 


=  5.25  *  12  *0.2836  /  2  [1 1]T 


f2  =  p  A  L  /  2  [1  1]T  =  3.75  *  12  *0.2836  /  2  [1  1]T 


GLOBAL  MATRICES 


K  =  3xl07  / 12  T  5.25  -5.25  0  1 

I  -5.25  9.00  -3.75  I 
L0  -3.75  3.75J 

F  =  [8.93  15.31+100  6.38]T 

BOUNDARY  CONDITIONS 

Qi  =  0 


MODIFIED  MATRICES 


K’=  3xl07/  12  T  5.25  -5.25  0  1 

I  -5.25  9.00  -3.75  I 
L0  -3.75  3.75J 

=  3xl07/12r  9  -3.75  1 

L-3.75  3.75J 

F’=  [8.93  (15.31+100-K21*a1)(6.38-K31*a1)]T 
=  [115.31  6.38]T 


Q’  =  [Q2  Q3]T 


FINAL  EQUATIONS 


K’Q’  =  F’ 

Solving  we  get: 

Q’  =  [0.9272  0.9953]  x  105” 

Q  =  [  0  0.9272  0.9953]  x  lO  5” 

STRESSES: 

a1  =  EB*q  =  3xl07  *  1/12  [1  -1][0  0.9272xl0  5]T 
=  23.18  psi 

a2  =  EB2q  =  3xl07*l/12  [1  -1]  [0.9273  0.9953]T*10 
=  1.70  psi 


